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ABSTRACT.  It  is  shown  tl  at  the  early  part  of  the  trajectory  of  ary  store  or  escape  capsule 
released  from  a  parent  aircraft  is  governed  principally  by  two  effects.  The  first  is  the  relative  velocity 
between  the  store  and  the  aircraft  at  the  instant  of  release.  The  second  is  the  acceleration  acting  on  the 
store  at  the  same  instant.  This  acceleration  is  due  not  only  to  gravity,  but  also  to  the  aerodynamic  forces 
and  moments  acting  on  the  store.  The  (normalized)  relative  velocity  and  (normalized)  acceleration,  as 
the  ordinate  and  ab  rissa,  define  a  planar  coordinate  system.  A  boundary  delineating  safe  from  unsafe 
separation  characteristics  can  be  drawn  on  this  plane.  Thus,  if  the  proper  data  are  available,  safe 
separation  conditions  can  be  predicted  in  advance.  The  available  data  from  flight  tests  strongly  support 
the  predictive  aspects  of  this  diagram. 

In  principle,  since  the  diagram  is  based  on  velocity  and  airloads,  no  additional  information  is 
needed  to  define  safe-separation  characteristics.  Moreover,  the  results  from  this  study  suggest  that 
airload  data  for  the  store  on  a  rack  car  be  directly  applied  to  give  valuable  results  about  safe-store 
separation. 
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' It  ia  shown  that  the  early  part  of  the  trajectory  of  any  store  or  escspe 
capsule  released  from  a  parent  aircraft  is  governed  principally  by  two  effects. 

The  first  is  the  relative  velocity  between  the  store  and  the  aircraft  at  the  instant 
of  release.  The  second  is  the  acceleration  acting  on  the  store  at  the  same  Instant. 
This  acceleration  is  due  not  only  to  gravity,  but  also  to  the  aerodynamic  forces 
and  momenta  acting  on  the  store.  The  (normalized)  relative  velocity  and  (normalized) 
acceleration,  as  the  ordinate  and  abcissa,  define  a  planar  coordinate  system.  A 
boundary  delineating  safe  from  unsafe  separation  characteristics  can  be  drawn  on 
this  plane.  Thus,  if  the  proper  data  are  available,  safe  separation  conditions  can 
be  predicted  in  advance.  The  available  data  from  flight  tests  strongly  support  the 
predictive  aspects  of  this  diagram. 

In  principle,  since  the  diagram  is  based  on  velocity  and  airloads,  no  additional 
information  ia  needed  to  define  safe-separation  characteristics.  Moreover,  the 
results  from  this  study  suggest  that  airload  data  for  the  store  on  a  rack  can  be 
directly  applied  to  give  valuable  results  about  safe-store  separation.  / 
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PREFACE 


In  recent  decades,  developments  in  store /carriage /launch  systems 
have  increased  tenfold  the  weapon-loading  combinations  available  to  a 
single  aircraft.  It  has  thus  not  been  possible  to  certify  the  launch 
and  separation  adequacy  of  all  possible  weapon  loadings  by  flight  test 
or  laboratory  investigation,  and  operational  problems  have  arisen  as  a 
result,  principally  with  free-fall  stores.  These  problems  have  prolif¬ 
erated  commensurately  with  the  growing  variety  of  free-fall  hardware, 
and  have  paced  in  severity  the  steady  increase  in  jet  aircraft 
performance . 

The  present  work  develops  a  general  separation  criteria  which  may 
first  find  useful  application  as  a  screening  device  for  quickly  sep¬ 
arating  those  cases  which,  by  the  criteria,  clearly  are  safe  from  those 
which  are  clearly  unsafe,  thus  quickly  identifying  that  vastly  reduced 
number  of  cases  which  are  marginal,  requiring  closer  scrutiny.  In  this 
way,  boundaries  of  release  envelopes  may  be  established  for  individual 
aircraft  with  greater  facility  and  confidence. 

In  Section  2,  the  criteria  are  developed  and  compared  with  the 
results  of  others.  A  detailed  mathematical  derivation  of  the  needed 
trajectory  information  is  given  in  Section  3,  which  can  be  skipped  with¬ 
out  loss  of  continuity,  if  desired.  A  summary  is  provided  at  the  end 
of  Section  2  that  includes  ’the  definition  of  the  several  parameters, 
and  provides  a  step-by-step  outline  for  the  application  of  the  proposed 
criteria. 
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LIST  OF  SYMBOLS 


¥ 


A  vertical  fall  parameters,  see  p.  14  (Eq.  24) 

Cj)  store  drag  coefficient 


store  axial  force  coefficient  in  aircraft  coordinates  and 
free  stream  conditions  at  time  t  ■  0  (+  forward) 


store  side  force  coefficient,  in  aircraft  coordinate  at 
time  t  *  0  (+  starboard) 


store  vertical  force  coefficient,  in  aircraft  coordinate 
at  time  t  *  0  (+  downward) 


store  lift  coefficient 

Cm'  store  pitching  moment  coefficient  about  aircraft  y  axis  at 

time  t  -  0  (+  nose-up) 


store  yawing  moment  coefficient  about  aircraft  z  axis  at 
time  t  »  0  (+  nose  starboard) 


6CX' — 6Cn'  change  in  coefficient  due  to  aerodynamic  interference. 

0  °  Note  Cx^  +  6CX^  is,  for  example,  the  coefficient  of  axial 

force  acting  at  the  instant  of  release.  The  same  applies 
for  the  other  coefficients. 


d 


2rmax  -  store  reference  length  for  moment  coefficient 


dr  relative  position  of  store  in  reference  system 


force 


) 


9 

*x»*y>*z 
kv ,  kz 


A  A  A 

i.J.k 


rack  impulse,  see  p.  2  (Eq.  3) 
applied  force 

acceleration  due  to  gravity 

moments  of  inertia  about  store  x,y,z  axis 

o 

store  radius  of  gyration  about  y,z  axis  Iy  -  mk^,  etc. 
unit  vectors  in  store  coordinates 


* 
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length  from  center  of  mass  to  3tore  nose 
&T  length  from  center  of  mass  to  store  tail 

m  store  mass 

Mq  slope  of  store  pitching  moment  coefficient 

[R]  transformation  matrix  from  store  coordinate  system  to 

reference  coordinate  system 

r*  distance  from  origin  of  reference  coordinate  system  to 

aircraft  center  of  mass 

r(x)  distribution  of  store  radius  in  store  coordinate 

rmax  maximum  value  of  store  radius 

O 

3  aerodynamic  reference  area  S  =  ffr 

t  time  seconds 

U(0)  store  velocity  in  x'  direction  at  t  =  0 

U'  store  velocity  in  reference  system  (primed  coord)  with 

respect  to  earth  (Section  2) 

-X 

V  reference  velocity  with  respect  to  earth  (Section  3) 

•A 

VA/C  aircraft  velocity  with  respect  to  earth  (Section  3) 

w  vertical  velocity 

W  store  weight 

x,y,z  store  coordinates,  origin  at  center  of  mass 

x',y',»'  reference  coordinates,  usually  origin  is  on  rack  or  pylon 


6 y'  initial  lateral  distance,  between  store  and  nearest  strike 

point 

<$z'(x,t)  distance  between  store  and  aircraft  in  reference  coordinates 
zc  displacement  due  to  constant  force 

6x,<5y,  etc.  displacement  in  x,y  direction 
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A3(x) 


A2(x) 


A3(x) 


A4(x) 


A5(x) 


Vx) 


relative  vertical  velocity  of  store  at  position  x  with 
respect  to  aircraft,  in  reference  coordinate  system  at 
t  •  0  (Section  2) 

relative  vertical  acceleration  of  store  at  position  x  with 
respect  to  aircraft,  in  reference  coordinate  system  at 
t  *  0  (Section  2) 

relative  lateral  velocity  of  store  at  position  x  with 
respect  to  aircraft,  in  reference  coordinate  system  at 
t  *  0  (Section  2) 

relative  lateral  acceleration  of  store  at  position  x  with 
respect  to  aircraft,  in  reference  coordinate  system  at 
t  *  0  (Section  2) 

relative  velocity  of  store  at  position  x  and  aircraft  in 
x  direction  in  reference  coordinate  system  at  t  *  0 
(Section  2) 

relative  acceleration  of  store  at  position  x  and  aircraft 
in  x  direction  in  reference  coordinate  system  at  t  **  0 
(Section  2) 


£  an  arbitrarily  small  parameter 

^'o* 0o*^o  store  Euler  angles  at  t  *  0 

*',0;,*'  aircraft  Euler  angles  at  t  -  0 

(ip*  0 * (p j  store  to  earth  coordinate  transformation 

p  dir  density  frr  from  aircraft  at  store  altitude 

AtUy,Awz  relative  angular  velocity  imparted  at  ejection  by  rack, 

t  *  0 


Note: 


angular  velocity  of  store  in  store  coordinate  system 

angular  velocity  of  aircraft  in  reference  system 

All  aerodynamic  force  coefficients  based  upon  free  stream 
dynamic  pressure,  q  and  S. 

All  aerodynamic  moment  coefficients  based  upon  free  stream 
dynamic  pressure,  q  and  Sd. 
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Section  1.  INTRODUCTION 


The  goal  of  the  study  reported  here  is  the  development  of  general 
criteria  to  enable  the  designer  or  operator  to  estimate  a  safe  store- 
separation  boundary  quickly  and  relatively  accurately.  Although  it  is 
likely  that  no  general  law  exists  stating  that  a  successful  aircraft  will 
be  put  to  an  increasingly  wide  set  of  tasks,  there  is  strong  empirical 
evidence  to  suggest  such  is  the  case.  Thus  an  aircraft  may  become  an 
all-purpose  carrier,  having  been  retrofitted  with  a  variety  of  multi¬ 
purpose  racks  and  pylons  for  carrying  an  astonishing  variety  of  shapes 
and  densities  of  external  stores.  In  case  of  emergency  the  pilot  must  be 
able  to  jettison  safely  these  stores,  pylons,  and  racks  in  a  great  variety 
of  combinations,  without  causing  damage  to  the  aircraft  from  the  store 
striking  the  aircraft.  From  an  operational  viewpoint,  the  safe-jettison 
envelope  should  be  larger  than  the  aircraft  operational  envelope.  Super¬ 
ficially,  the  safe-jettison  requirement  would  seem  to  be  the  most  strin¬ 
gent  of  all  the  store-separation  requirements.  Actually,  this  is  not 
necessarily  true.  A  low-drag,  low-density  store,  i.e.,  a  droppable  fuel 
tank,  or  an  empty  gun  pod  suspended  near  a  general-purpose  rack,  may 
represent  a  more  difficult  store-separation  problem. 

These  problems  of  store  separation  have  received  considerable 
attention  recently.  Black  (Ref.  1)  has  shown  that  in  some  cases  wind 
tunnel  tests  are  extremely  useful  in  predicting  full-scale  separation 
characteristics.  Earlier  work  has  been  summarized  by  Reed,  et  al 
(Ref.  2).  The  question  of  general  scaling  has  been  widely  discussed, 
and  a  summary  of  scaling  for  free  drop  tests,  the  grid  method,  and 
trajectory  simulation  is  discussed  in  Ref.  3  and  4,  Application  of 
these  results  to  A  specific  problem  is  contained  in  Ref.  5,  as  is  a 
suggestion  for  criteria.  D.  A.  Jones  (Ref.  6)  has  presented  a  straight¬ 
forward  analysis  of  possible  criteria.  These  criteria,  and  Unger's 
suggestions  (Ref.  7),  are  discussed  herein  in  Section  2. 

In  establishing  useful  design  criteria,  it  is  necessary  to  be  con¬ 
servative  to  allow  for  variations  in  geometry,  performance,  and  opera¬ 
tional  conditions,  while  still  defining  a  safe  separation.  Necessarily, 
several  stringent  assumptions  are  made,  but  the  evidence  presented 
herein  supports  their  validity.  These  include  assumed  aerodynamic  loads 
arid  rack  ejection  characteristics.  The  former  may  either  be  measured  or 
estimated  analytically.  These  computations  must  include  ways  or  account¬ 
ing  for  not  only  the  far  field,  slowly  varying  aerodynamic  interference 
between  the  store  and  the  aircraft,  but  more  importantly  the  near  field, 
rapidly  varying  aerodynamic  interference.  The  far  field  problem  resolves 
to  that  of  determining  the  aerodynamic  characteristics  of  the  store  in 
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a  nonuniform  field.  Results  from  a  comparison  between  calculations  and 
wind  tunnel  tests  suggest  that  slender  body  theory,  scaled  to  the  free- 
air-uniform  field  experimental  results,  and  expressed  as  follows,  is 
suitable 


Lift 


1  Cl<*  MEAS  '  Sre f/*  ,  dS 

"  2  pu  "s(£N)  -  s Ur)  j  a'x,y,z)  dx 


dx 


w  _  1  „,,2  ^  MEAS  SreF^REF  I" 

Moment  ■  —  pU^ - 1  < 


Volume  •  (ko 


o(x,y  ,2)  (x,- 


.  dS 


(1) 


(2) 


The  notation  a(x,y,z)  implies  that  the  flow  curvature  is  spatially 
dependent.  The  remaining  notation  is  standard. 

The  near  field  interference  has  yet  to  be  treated  successfully  in 
an  analytical  way.  In  the  general  case,  this  interference  is  due  not 
only  to  potential  flow  interaction  between  stores  and  wing,  stores  and 
pylon,  or  stores  themselves,  but  also  to  a  viscous  flow  separation  and 
sometimes  shock  wave  that  acts  to  repel  the  stores.  Usually  near  inter¬ 
ference  is  of  an  attractive  nature;  its  decay  rate  is  strongly  shape- 
dependent.  In  the  case  of  a  sphere  the  near  field  interference,  pro¬ 
portional  to  (rmax/dz’)^,  is  down  to  0.6  of  its  initial  value  in  the 
distance  of  one  diameter.  These  near  interference  effects  usually  act 
over  a  short  time  period,  say  0.2  second  or  less.  Their  time  scale  is 
the  same  as  the  time  scale  of  the  rack  ejector  mechanism. 

A  typical  characteristic  of  ejector  racks  is  shown  in  Fig.  1  (from 
Fig.  13  of  Ref.  3).  While  this  characteristic  is  not  that  of  a  true 
impulse,  which  implies  a  jump  in  velocity  while  a  jump  in  force  implies 
a  jump  in  acceleration,  it  can  be  so  treated.  Note  this  force  terminates 
at  about  0.06  to  0.10  second.  It  will  be  assumed,  in  the  analysis,  that 
the  actual  impulse  will  be  applied  at  zero  time  and  zero  relative  motion. 
Some  care  is  needed  at  this  point  because  a  flexible  rack  may  deflect 
under  the  impulse  rather  than  accelerate  the  store.  For  this  calculation 
a  force  distance  characteristic  is  needed,  but  with  a  reasonable  assump¬ 
tion  the  force  distance  curve  may  be  converted  into  a  force  time  curve. 


The  structure  of  these  equations  is  based  upon  linear  theory  (cf 
Nielsen,  J.  N.,  Missile  Aerodynamics,  McGraw  Hill  1960,  pp.  114-118, 
182-201);  experimental  support  is  based  upon  the  results  of  Beane  and 
Durgin  (ASD-TDR-61-295  part  2). 
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The  near  field  interference,  as  can  be  shown,  aleo  has  an  effect 
that  acts  like  the  force  distance  effects  discussed  above.  This  ^gain 
differs  from  a  true  impulse.  However,  we  dll  regard  the  near  field 
effects  as  a  constant  force  and  moment  acting  from  the  release  point. 

The  primary  problem  is  to  estimate  the  size  of  near  field  effects  with 
sufficient  accuracy. 

The  criteria  that  will  ba  developed  may  be  shown  on  a  plane  coord¬ 
inate  system  whose  ordinate  is  the  relative  velocity  of  any  selected 
point  on  the  surface  of  the  store  (i.e.,  at  x)  with  respect  to  the  rack, 
in  a  rack  coordinate  system.  The  abcissa  is  the  acceleration  of  the 
same  point,  i.e.,  the  gravity  component  plus  that  caused  by  airloads 
acting  at  the  instant  of  release.  It  turns  out  that  such  planes  are 
needed  for  both  vertical  and  lateral  motion  (and  axial  motion  for  pilot 
escape  capsules  and  other  special  cases) .  It  is  possible  to  draw  a 
line  or  lines  on  these  planes  that  separate  safe  from  unsafe  store 
spearation.  The  validity  of  the  proposed  criteria  has  been  tested 
against  flight  test  data. 

From  a  pragmatic  point  of  view,  a  miss  of  1  foot  or  one  radius  may 
be  considered  the  limit.  Rather  than  specify  the  details  of  the  rack 
and  pylon,  we  shall  pass  a  plane  through  the  lowest  point  of  the  rack 
parallel  to  the  aircraft  lateral  plane  and  define  this  plane  as  the 
critical  plane,  Then,  any  store  that  fails  to  move  one  radius  away  and 
normal  to  this  plane  in  0.25  second  will  be  assumed  to  be  not  safe^ 

(Ref.  5).  This  latter  requirement  will  be  expressed  more  rigorously  in 
terms  of  characteristic  times  and  distances,  but  this  definition  seems 
to  be  adequate  for  the  purpose  of  this  report. 


The  choice  of  one  radius  in  1/4  second  is  the  author's.  The 
basic  idea  is  David  L.  Schoch's, 
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Section  2.  SAFE-SEPARATION  CRITERIA 


The  criteria  developed  for  safe  store  separation  will  be  discussed 
in  detail  in  this  section,  which  also  contains  available  data  in  support 
of  the  usefulness  of  the  criteria.  The  detailed  mathematics  supporting 
the  choice  of  variables  used  to  define  the  safe-separation  criteria  (SSC) 
is  presented  in  Section  3. 

The  set  of  variables  used  to  define  the  SSC  is  based  upon  one  pri¬ 
mary  assumption:  the  trajectory  near  the  aircraft  is  governed  by  the 
forces  and  moments  acting  on  the  store  at  the  instant  of  launch  and  they 
continue  to  act  at  the  same  magnitude.  This  assumption  allows  simple 
computation  of  a  trajectory  in  terms  of  dimensionless  variables.  Once 
the  proper  variables  arc  selected,  a  limiting  boundary  in  a  plane  defined 
by  these,  variables  can  be  drawn  either  by  theory,  test  data,  or  some 
combination  of  both.  While  the  variables  used  in  the  SSC  presented  here 
are  based  upon  a  constant  force  and  moment  approximation  to  the  actual 
data. 


To  validate  the  constant  force  assumption,  one  can  compare  actual 
trajectories  with  trajectories  based  upon  constant  forces  and  moments. 
The  time  dependence  of  vertical  distance  as  a  solution  of  a  second- 
order  differential  equation  for  z  is 


z(t)  =  ^ Jfl  F(T)dTdt 


(3) 


in  response  to  an  unsteady  force  F(t) .  Suppose  F(t)  can  be  represented 
by  a  Taylor  series 


F(t) 


F(0)  +  l 
n=l 


dnF 

dTn 


^n 

n! 

t  =  0 


(4) 


Then 


z(t) 


1 

F(0;t2 

00 

+  I 

dnF 

tn+2 

m 

2! 

n=l 

dtn 

(n+2)! 

t  =  0 


w(0)t  +  z  (0) 


(5) 
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Clearly  Che  effects  of  variation  of  the  force  upon  the  displacement 
always  appear  two  powers  of  t  higher  than  in  the  series  for  the  force. 
Thus  when  the  value  of  z  for  the  constant  force  approximation  (zc)  is 
arbitrarily  close  to  the  actual  z,  (|z-zcj<  e) ,  there  is  corresponding  to 
each  c  a  sufficiently  small  value  of  t.  The  practical  problem  is  whether 
or  not  this  small  value  of  t  is  large  enough  in  practical  cases.  Detail 
examination  chows  the  initial  force  is  the  difference  of  gravity  and  an 
aerodynamic  force.  Hence  the  constant  force  assumption  produces  trajec¬ 
tories  that  lie  above  (smaller  z)  the  actual  trajectories.  The  actual 
variation  of  forces  and  moments  in  two  cases  are  shown  in  Table  1.  It 
can  be  seen  that  the  variation  of  the  forces  is  as  high  as  a  factor  of 
30  in  the  first  0.2  second.  Nevertheless,  the  vertical  trajectories  are 
in  fairlv  close  agreement  (Fig.  2).  The  lateral  trajectories  seem  to  be 
more  sensitive  since  the  lateral  forces  are  small  to  begin  with. 


TABLE  1.  Variation  of  Forces and  Moments. 


Time , 

fn. 

lb 

Fy 

,  lb 

M,  lb -ft 

N,  lb-ft 

sec 

66 

67 

66 

67 

66 

67 

66 

67 

0 

913 

}  5 .4 

-25.7 

-43.0 

-30.6 

-38.5 

-99.6 

-1.34 

0.05 

33.7 

31.7 

-33.8 

-53.9 

-19.6 

-20.6 

-86.7 

-1.10 

0.10 

12.6 

-152.0 

-61.5 

-1.25 

-927.8 

171.3 

-58.2 

-72.3 

0.15 

-579 

53.6 

-93.4 

1.94 

-785.3 

23.1 

-30.8 

-18.2 

0.20 

-35.9 

-1.0 

-1.25 

-2.6 

-633.0 

30.4  j 

-724.0 

37.81 

It  is  remarkable  that  a  position  error  of  only  33%  is  introduced  by 
the  constant  force  approximation  while  (for  example)  the  force  varies  from 
913  to  -36  lb.  For  this  reason  the  parameters  were  selected  on  the  basis 
of  the  constant  force  assumption.  Figure  2  shows  that  the  constant  force 
assumption  generally  underestimates  the  z  displacement,  which  is  conserva¬ 
tive,  and  also  underestimates  the  y  displacement,  which  is  not  conservative; 
but  a  method  of  allowing  for  lateral  errors  on  the  safe-separation  envelope 
will  be  given  below.  The  final  SSC  is  felt  to  be  conservative. 


INITIAL  DEVELOPMENT  OF  THE  SAFE-SEPARATION  CRITERIA 

The  SSC  is  developed  from  considerations  of  the  distance  between 
some  point  on  the  store  and  a  reference  point  in  the  critical  plane 
described  in  the  Introduction.  The  vertical  distance  is  emphasized  here 
because  with  a  suitable  time  constraint  it  can  be  made  to  be  a  sufficient 
condition.  Schoch  (in  a  private  communication)  argues  that  the  store 
should  quickly  fall  a  critical  distance  equal  to  the  maximum  diameter  of 
the  store.  In  that  case  a  store  which  moved  sideways  to  the  wing  tip 
without  striking  the  aircraft  and  then  fell  one  diameter  would  be  considered 
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FIG.  2.  Comparison  of  Constant  Force  Trajectories  With 
Actual  Trajectories. 
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unsafe^.  The  initial  vertical  distance,  6z'(0),  is  the  distance  between 
the  Initial  contact  point  and  the  point  on  the  store  in  question  in  a 
coordinate  system  fixed  to  the  rack;  2'  is  thus  an  aircraft  fixed  coordi¬ 
nate,  positive  downward.  Generally  6z'  is  a  function  of  x  and  t,  i.e., 
(from  Section  3,  Eq.  93) 


6zf (x,t) 


<5z '  (0,0)  +  Ax(x)t  +  A2(x)2“  + 


(6) 


BASIS  OF  THE  CRITERIA 

To  develop  the  SSC  we  want  to  find  out  the  circumstances  for  which 
Eq.  6  has  no  real  positive  roots,  because  a  positive  root  implies  that 
the  store  hits  the  aircraft.  Setting 

6z'(x,t)  >  6z' (0,0)  (7) 

gives  for  the  quadratic  (constant  force)  approximation 

t(~  t  +  AX)>  0  (8) 

Thus  the  roots  are  t  =  0  and  t  =  -2AjyA2.  If  A^  and  A2  each  are  greater 
than  zero,  there  will  never  be  a  positive  root,  so  the  separation  is 
safe.  Ai  and  A2  are  defined  to  be  the  initial  velocity  and  initial 
acceleration  just  at  release  of  the  store  in  the  z'  direction  (note  the 
case  A^<0  and  A2<0  has  been  excluded  for  stores  suspended  below  the 
aircraft)  .  Hence  we  can  say  that  ApO  and  A2>0  is  a  sufficient  condition 
in  the  z'  direction  for  SSC.  Ai  and  A2  have  dimensions.  Usually  non- 
dimensional  parameters  are  preferable,  sc  we  use  the  characteristic  velocity 


yj  2^rmax 


and  the  acceleration  due  to  gravity,  g.  Defining  a  plane  coordinate 
system  with  the  ordinate 

Al  j \j  2g,rmax 

and  the  abcissa  A2/gf,  we  can  say  that  a  store  whose  initial  velocity  and 
initial  acceleration  fall  in  or  on  the  axes  of  the  first  quadrant  of  this 


The  lateral  case  will  be  discussed  later. 
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plant  has  satisfied  the  safe-separation  criteria.  This  is  indicated  in 
the  sketch  below.  The  line  Ai  *  C  is  included  in  the  safe  region  on 
the  right  half  plane. 


l  •r  m«x~ 


8CMOCWS  LIMIT  DEFINED 
•Y  FALLING  TOO  SLOWLY 

/V 


// s 


>  ; 


//>. 


SAFE  REGION 
(VERTICAL  DR  OF) 


^2 


I 


The  useful  region  of  the  &i“A2  plane  (the  first  quadrant)  is  further 
restricted  by  time  arbitrary  requirement  that  the  store  fall  one  maximum 
radius  in  one-qparter  second.  This  limit  line  is  similar  to  the  criteria 
used  by  Schoch?  Thjis  new  limit  is  derived  from  the  distance  relation, 
which  may  be  W£^r£en 


(9) 


Setting  6z'  =  rmax  and  normalizing  gives 


A1 

7^ 


1  6z 
yj~2  rmax 


'  .  /rmax 

ax  V  t^g 


j,- 

2'J~2.  rroa*  \  9  J 


(10) 


The  term  V  rmax/t2 g  times  the  ^2  is  the  ratio  of  the  time  for  the  store 
to  fall  its  own  radius  at  one  g  to  the  critical  time.  This  relation 
between  A^  and  A2  defines  a  diagonal  line  of  negative  slope  that  excludes 
the  region  near  the  origin.  This  line's  location  depends  upon  the  radius 
(i.e.,  a  store)  and  a  critical  time.  For  a  12-inch  store  and  a  quarter 
second,  for  example 
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^1  /  ^2 

—  £=  =:  0.35  -  0.71  (  — 

The  broken  line  shows  this  boundary  on  the  sketch  above.  It  Is  seen  that 
the  restriction  is  not  large.  Indeed,  from  the  example  calculated,  the 
excluded  region  corresponds  to  a  small  triangle  bounded  by  a  velocity  of 
about  35%  of  that  attained  in  free  fall  and  by  an  acceleration  of  0.49  g. 
This  acceleration  is  equivalent  to  an  airframe  pitch  angle  of  greater  than 
60  degrees,  a  steep  climb  or  dive.  Note  further  that  the  slope  of  this 
line  is  proportional  to  the  critical  time  and  the  intercept  is  propor¬ 
tional  to  the  reciprocal  of  the  critical  time.  Hence,  the  shorter  this 
time,  the  larger  the  region  of  the  first  quadrant  that  is  excluded  from 
the  region  of  operation.  Note  that,  as  the  location  of  the  characteristics 
of  the  store  on  this  plane  are  farther  from  the  boundaries,  the  cleaner 
the  separation  characteristic. 

Briefly,  it  seems  possible  to  define  a  boundary  in  the  velocity- 
acceleration  plane  that  delineates  those  stores  that  will  separate 
safely  from  those  that  will  not.  Ultimately,  the  best  boundaries  should 
be  determined  by  experimental  results . 

Figure  3  shows  typical  plots  on  this  vertical  SSC  plane.  The  solid 
line  is  taken  from  computations  by  Daniel  Jones  of  the  Naval  Weapons 
Laboratory  (Ref.  6).  This  line  lies  to  the  left  of  the  ordinate  because 
the  negative  acceleration  effects  considered  become  smaller  as  the  store 
moves  away  from  the  aircraft.  Hence,  high  ejection  velocities  can  allow 
safe  clearance  for  the  store  in  the  left  half  plane  unless  they  induce 
strong  restoring  aerodynamic  forces  that  cause  the  store  to  rise  back 
to  the  aircraft. 

Before  discussing  Ai  and  A2  in  detail  let  us  consider  the  lateral 
case.  Here  we  have  two  alternatives.  The  first,  simplest,  and  most 
conservative  alternative  is  to  simply  let  the  store  move  almost  the 
distance  6y'(0)  to  the  nearest  store.  Again,  from  Section  3  we  find 
this  implies 
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In  this  case,  depending  upon  the  direction  of  the  nearest  store,  either  f 

the  first  or  third  quadrants  are  acceptable,  as  sketched  below. 


Note  the  first  quadrant  applies  if  the  store  must  move  to  the  left 
(6 y*(t)  <  0)  to  strike  something  solid,  while  the  third  quadrant  applies 
to  the  case  where  movement  to  the  right  is  dangerous.  However,  there  is 
a  better  alternative  to  define  a  boundary  in  this  plane.  This  boundary 
comes  about  because  the  lateral  motion  of  a  store  may  be  limited.  Hence, 
we  will  require  a  store  that  moves  a  distance  6y'  (the  distance  to  the 
object  it  may  strike)  must  also  fall  a  distance  6z'  in  tne  same  time. 

6z'  is  the  distance  the  store  must  fall  to  pass  under  the  object  at  5y'. 
Arguing  as  above  (Eq,  9  and  10) 


Now,  t,  the  time  to  fall  the  distance  5z' ,  i3 


t 


(14) 
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Substituting  and  simplifying, 

A3  „  6y'  A  \ 

J  2ar  2rmax 
V  ‘9rmax 


where 


q  '/23Tmax  /  A1 
A  =  —a— - — ;  ===  ——  ■  ■  ■  — - 

\fl  +  §*—  .  .  2grmax  _  , 

V  rMX  g  A2  1 


(15) 


(16) 


Note  if  Ajl  ■+  0,A  ->  1. 

This  line  in  the  A3,  A4  plane  has  a  negative  slope  and  limits  the  useful 
area.  Thus  the  safe-separation  region  in  the  lateral  plane  is  determined 
by  two  slanted  lines.  The  intercept  in  the  A3  axis  is  determined  by 
6y'(0).  The  larger  the  distance  Sy'(O),  the  bigger  the  useful  area.  For 
many  practical  calculations,  62 '  can  be  taken  to  be  one  diameter  or  one 
fin  span.  Note  this  boundary  is  fixed  by  a  particular  store-aircraft 
situation.  This  latter  alternative  seems  to  be  the  better. 

So  far  the  fore  and  aft  motion  of  the  store  has  not  been  considered. 
In  every  case  studied,  a  store  whose  fore  or  aft  motion  causes  trouble 
was  in  trouble  from  the  vertical  or  lateral  criteria.  However,  if  one 
writes,  for  completeness,  an  x  displacement  equation 


fix' 


A5t  +  A6 


£ 

2 


(17) 


with  A5  usually  zero  and  A^  <  0.  The  nondimensional  relation  insuring 
that  the  store  falls  farther  than  it  moves  aft,  then,  is 


fix' 

< 

i. 

rmax 

2 9  A2 

(18) 
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(l/2)pU2(0)SCD 


This  term  may  be  large  if  the  store  rou  tes  90  deg,  but  as  indicated 
above,  usually  a  store  that  rotates  90  deg  this  close  to  the  aircraft  is 
already  suspect.  Alternately,  a  line  can  be  drawn  on  the  A^,  A2  plane  to 
insure  that  this  limit  be  met.  Using  earlier  results 


6z' (t)-6z|(0) 


and  from  the  x  relation 


Substituting  and  nondimensionalizing  gives  the  relation 


!  A,r 
6  max 


^rmax 


,r 

6  max 


In  the  relation,  Ag  and  6x*  are  negative.  This  line  represents  a  second 
lower  bound  and  may  be  compared  with  the  earlier  condition,  where  the 
intercept  is 


as  compared  with 


^6  rmax 


2rmax  V  9  fix' 


Similarly,  the  slopes  are 


yw. 

V  Ac,rmax 


14 


'«!  #*rt««**»*'»*'»Wh**'  w  r  - 


NWC  TP  4995 


Since  <Sx*  is  usually  large,  Cq  is  small  and  6z*  is  small,  the  aft  distance 
criteria  lies  below  the  time-to-fall  criteria  as  long  as 

6x*  t£ 

*6  2 

Hence  this  criteria  is  not  usually  needed,  although  it  may  apply  to 
items  like  safe/arm  pins  and  arming  wires. 


APPLICATION  TO  PILOT  ESCAPE  CAPSULES 

Application  of  this  result  to  the  pilot  escape  capsule  requires 
some  minor  changes,  mostly  in  definitions.  Since  the  capsule  leaves  from 
inside  the  aircraft,  A^,  A2,  Ag  will  be  defined  at  the  instant  the  capsule 
just  clears  the  fuselage.  This  instant  will  be  taken  to  be  the  origin 
for  time,  i.e.,  t  =  0  at  this  instant.  For  an  upward  ejection  Aj  <  0, 

A2  >  0  and  surely  Ag  <  0.  If  the  highest  point  on  the  aircraft  has 
coordinate  -x^'  -z^1  then  the  safe  region  in  the  Aj,  A2  plane  will  be 
mostly  in  the  third  quadrant,  below  the  line 


Here  h  is  the  maximum  height  of  the  capsule.  Thus  the  safe  region  for 
the  upward-ejected  pilot  escape  capsule  is  below  and  to  the  left  of  the 
line  as  shown  below. 


NWC  TP  4995 


Note  that  even  if  A2  >  0,  there  is  a  domain  in  the  fourth  quadrant  that 
allows  safe  operation.  Previous  capsules  probably  operated  in  this 
domain. 


RELATION  TO  EARLIER  RESULTS 

All  the  safe-separation  criteria  discussed  heretofore  are  based 
upon  the  idea  of  displacement,  either  in  the  complete  form  given  in 
Section  3  or  in  one  or  another  simplification. ^  Jones,  Schoch,  and 
Unger  adopted  the  simplified  form.  (Schoch  has  used  a  better  criteria 
by  using  actual  stores  and  actual  aircraft  shapes.)  In  general,  the 
simplified  spatial  relationship 

z'(t)  -  «(t)  >  zAC(t)  (24) 

is  employed  where  l  >  0  for  the  nose  and  l  <  0  for  the  tail.  A  slightly 
more  accurate  form 


z'(t)  -  l  sin  8  -  t(Z)  cos  0  >  zA(,(t) 


(25) 


simply  allows  for  the  possibility  that  there  is  a  store  shape  change.  In 
each  case  that  a  simple  result  is  sought,  as  opposed  to  a  detailed  numer¬ 
ical  result,  the  aircraft  is  represented  as  a  plane.  Schoch  adopted  the 
reference  plane  below  the  store  so  that  6z'(0),  the  initial  store  refer¬ 
ence  plane  distance,  is  negative. 


The  treatment  of  this  basic  formula  (Eq.  25)  varies  from  investi¬ 
gator  to  investigator.  Jones  calculates  the  displacement  at  time 

t  =  tt  >/  Iy/  -  although,  in  the  same  vein,  he  could  have  selected 

t  =  2  v  ly/  ~  Mot*  Jones  then  used  the  constant  force  equation  evaluated 


at  the  half  period,  or  the  condition  of  finite  downward  velocity  at  the 
edge  of  the  lift  interference  field. 


Unger  used  a  more  complete  solution  in  the  basic  condition  for  safe 
separation  and  simply  computed  trajectories,  isolating  those  that  struck 
the  aircraft.  He  also  expanded  his  solutions  to  the  4th  order  in  time 
and  substituted  into  the  equation  that  represents  the  separation  distance, 
as  well  as  writing  two  approximate  forms  involving  trignometric  force 


The  use  of  constant  force  and  moment  analysis  has  been  proposed 
also  by  G.V.R.  Rao  and  B.  Jaffee  "Study  of  Generalized  Safe-Launch 
Bounds  for  External  Stores,"  submitted  to  NSRDC  Aerodynamics  Laboratory, 
October  20,  1967,  on  Contract  N00600-67-C-0592 .  Rao  and  Jaffee  construct 
the  safe  limits  in  th«  P,  z  plane. 
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functions.  The  simplest  is  a  constant  force  approximation  (the  inter¬ 
ferences  do  not  vary) ,  thus 

-£Awy  C08  f-^Ji  A 

y  ypvijT  %  V  7 

+  wQt  +  F(0)'|—  >  zAC  (26) 


reduces  to  the  second-order  result  described  abov...  That  is,  we  have 
replaced  sin  x  by  x  and  cos  x  by  l-x^/2.  This  suggests  that  the  quad¬ 
ratic  approximation  indicates  the  store  will  rotate  too  fast  (cos  x  ■  0 
at  x  -  y/2  not  tt/2.)  The  implication  is  that  the  quadratic  approximation 
is  conservative,  because  the  store  falls  relatively  more  slowly  than  it 
rotates.  Hence,  any  store  that  is  safe  by  the  quadratic  or  constant 
force  criteria  is  likely  to  always  be  safe. 

Figure  3  shows  the  SSC  computed  from  the  constant  force  approxima¬ 
tion  and  a  typical  lower  boundary  due  to  distance-time  limits  and  the 
results  from  Unger  and  D.  A.  Jones.  The  most  striking  point  is  the 
difference  between  the  ordinate  as  a  boundary  and  a  boundary  derived 
from  Fig.  7  in  D.  A.  Jones*  report.*  Equally  striking  i9  the  effect 

*NOTE:  A  linear  approximation  gives  the  Jones  boundary  the  form,  for  a  stable  store. 


J  2?r 


max 


/2M(0) 


ir/-  ^(0)1 


1  T_  ^2 

2  tc  g 


Here  T  is  the  period  of  the  store  in  a  uniform  stream,  tc  is  the  critical  time.  This 
line  has  a  slope  that  is  much  steeper,  in  the  ratio  T/2tc,  than  Schoch’s  line  and  is 
the  boundary  for 


/ rmax  x 

/■M(O)  /  tt2v.  A“y 

y  1 2(7  j 

V  V9*  max 

y-Moi '  47  yi 

0.354 


yn l. 

rmax 
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Note  this  line  (because  of  the  ratio  M(0) / V -Ma)  is  not  fixed,  but  rather  depends 
upon  the  flight  conditions  such  as  q,  Mach  number  and  angle  of  attack.  The  latter  is 
dependent  upon  load  factor.  These  terms  are  known  because  they  are  contained  in  A2 , 
but  nevertheless  this  is  a  movable  boundary.  For  T  >>  tc  the  ratio 


^2(1  -  J-)x 


t  M(0) 
c 


9 


-  2.29 


n  V  g rraax  I 


^£?rmax 


xM(0) 

1 


For  T  «  tc,  the  limiting  value  of  A2 Ig  approaches  the  value  2rmax/t^g.  The  former 
case  is  more  likely  to  be  encountered. 


17 


NWC  TP  4995 


of  the  (one  radius  in  1/4  second)  limit  in  increasing  the  required  ejec¬ 
tion  velocity  in  the  region  of  negative  acceleration.  The  region  in  the 
second  quadrant  dominated  by  Jones'  results  are  those  cases  where  a  stable 
store  moves  down  and  then  up  to  hit  the  aircraft.  The  lack  of  symmetry 
in  the  vertical  plane  (A3,  A2)  when  compared  with  the  appearance  of  the 
A3,  A^  plane  (Fig.  3  and  4)  is  due  to  the  geometry  of  the  lateral  configur¬ 
ation  if  they  are  not  mounted  on  a  wing  tip  or  a  side  body  pylon. 

Returning  to  the  expressions  of  A3,  Aj,  A3,  and  A^  in  detail,  one 
finds  that  much  of  the  experience  gained  in  the  area  is  supported  by  this 
SSC.  There  are,  however,  several  features  that  are  worth  a  detailed 
discussion.  A3  and  A3  are  simple,  i.e.,  (from  Section  3) 

A3  ■  w  -  AtOy  x(cos  tj>0  sin  90  sin  0q  +  cos  0O  cos  8q) 

-  AWy  r(x)(cos  ifi0  cos  0O  sin  G0-  sin  0O  cos  8q)  (27) 

Thus,  A3  is  just  the  velocity  imparted  to  a  store  at  position  x  in  the 
coordinate  system  attached  to  the  rack,  w  is  the  vertical  ejection 
velocity  of  the  store  with  respect  to  the  aircraft  rack  reference  point, 
acting  at  time  t  *  0.  Similarly,  Aoiy  is  the  store  angular  velocity 
about  the  y*  axis  with  respect  to  the  aircraft  rack,  w  is  positive  if 
it  acts  downward  and  Au>y  is  positive  nose-up.  x  is  positive  for  stations 
ahead  of  the  center  of  mass  and  negative  for  stations  aft  the  center  of 
mass.  r(x)  the  radius  (or  fin)  distribution  is  always  positive.  By  the 
same  arguments, 

A3  -  v  +  Auz  [x  -  r(x)]  (28) 

From  Section  3,  the  acceleration  term  is,  for  straight  and  level 
flight, 


(X-  ^-r)-d-(CM;  +  5cm;)]  >  0  (29) 


In  this  relation,  Czq,  6CZg,  Cm^  and  6Cm^  are  evaluated  at  the  Mach  number, 
angle  of  attack  appropriate  to  the  load  factor,  and  angle  of  side  slip 
that  are  imposed  by  the  trajectory.  The  altitude  and  flight  speed  are 
characterized  by  the  dynamic  pressure,  and  the  aircraft  trajectory  appears 
in  cos  0q  cos  t>o. 


=  cos  0q  cos  <j>Q 


+  \  pu2<°>!  [ c*; +  5c*;  - 


18 


NWC  TP  4995 


The  following  conclusions  may  be  drawn  for  stores  below  the  aircraft: 

1.  The  effectiveness  of  gravity  is  reduced  for  any  non-straight 
and  level  flight  path. 

2.  If  the  term  inside  the  brackets  is  negative,  i.e»,  the  aerodyna¬ 
mic  loads  force  the  store  towards  the  aircraft,  a  lightweight  store  is 
more  likely  not  to  clear  the  aircraft.  This  situafton  is  aggravated  as 
(l/2)pU2(0)  becomes  larger. 

3.  If  the  term  inside  the  brackets  is  positive,  a  heavy  store  or 
a  low  dynamic  pressure  flight  is  more  sensitive  to  flight  path  angle, 
because  the  aerodynamic  forces  and  moments  are  less  effective. 

4.  The  constant  force  and  moment  approximation  suggests  that  the 
same  information  needed  to  evaluate  the  carriage  air  loads  on  a  store 
can  also  be  used  to  determine  the  safe-separation  criteria.  Thus  extra 
information  is  not  needed ,  but  the  needed  information  may  be  shared  by 
different  groups  of  the  same  design  team. 


PROPOSED  SSC 

The  SSC  is  developed  on  the  nondimensional  velocity-acceleration 
planes  by  applying  every  bit  of  data  the  author  could  find  for  which  all 
the  needed  parameters  were  available  This  data  is  listed  in  Table  2. 

The  values  of  Ai,  A2,  Aj,  and  A4  were  computed  for  both  the  nose  and  the 
tail,  and  the  results  plotted  together  with  the  limit  lines  for  the 
particular  store. 

Figures  5  and  6  show  the  location  of  nose  and  tail  points  on  the 
nondimensional  velocity-acceleration  plots  (the  run  numbers  are  contained 
in  the  circles).  The  data  is  for  conditions  listed  in  Table  2  of  Ref.  1. 
All  these  separations  were  clear  in  the  vertical  plane.  The  implication 
is  that  the  ordinate  is  far  too  restrictive,  and  that  even  the  1/4-second, 
one  radius  limit  is  too  severe.  The  nose  points  all  show  a  positive 
acceleration  but  almost  no  vertical  velocity  at  ejection,  whereas  the 
tails  come  down  because  of  their  velocity,  even  though  they  are  accelerated 
negatively  (upward) .  The  tails  of  runs  68  and  particularly  71  tend  to 
float,  i.e.,  it  takes  them  a  longer  time  to  clear  the  area.  In  the 
lateral  plane,  run  67  is  out  of  the  safe  domain — its  nose  moves  too 
far  to  the  left  (hits  the  fuel  tank).  Ir.  run  71  the  tail  brushes  the 
pylon.  The  analysis  says  it  is  safe  (A ^/g  'v  0.03)  but  between  "hanging 
up"  and  the  errors  in  the  measurement  (at  least  that  large)  the  clearance 


In  particular  the  author  is  indebted  to  Charles  Matthews  (AFATL/ATII 
Eglin  Air  Force  Base) . 
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TABLE  2.  Summary  of  Store  Conditions  at  Release 


Vertical  Separation  Conditions  for  Data  From  Ref 


DATA  FROM  AEDC  TR-69-130 
TABLE  JL  NUMBER  IS  RUN 
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is  uucertain.  Run  68,  wnich  i3  run  67  at  &  reduced  angle  of  attack 
(Men  -  0.7),  cleared  safely.  Figures  5  and  6  would  support  these  results 
if  Jones’  criteria  and  the  1/4-second  Unit  were  used. 

Figures  7  and  8^  show  the  characteristics  of  the  Mk  4  gun  pod  os  an 
A-6  pylon.  These  characteristics  are  based  upon  flight  test  trajectories 
and  computed  trajectories.  The  lateral  separation  was  clean  for  these 
cases.  On  test  533  the  tall  hit  the  pylon,  while  on  tests  553  and  556 
the  nor-e  was  within  inches  of  hitting  the  pylon.  The  plot  shows  the  nose 
on  553  and  556  actually  had  an  upward  velocity  at  release,  so  it  is  not 
surprising  that  it  hit.  The  tail  hit  on  test  533  is  also  Indicated, 
x'hese  flight  test  results  again  indicate  the  utility  of  the  acceleration- 
velocity  plots. 

Figure  9^  shows  the  vertical  separation  characteristics  of  the 
Sargent-Fletcher  fuel  tank  (empty)  from  the  0V-1D.  In  every  case  the 
trajectory  is  clear  and  the  plot  of  the  flight  test  results  shows  that 
jettison  is  safe.  Similarly  Fig.  10  shows  the  lateral  characteristics 
are  satisfactory.  Flight  57  is  most  Interesting,  because  here  the  air¬ 
craft  was  rolled  and  the  basic  values  of  A4/9  were  negative.  However, 
because  of  the  roll,  a  gravity  component  acted  to  cause  Che  store  to  fall 
in  the  outboard  direction.  Note  that  even  ?.  slight  negative  value  of 
t^/g  is  not  serious,  because  there  is  8  ft  of  distance  before  there  is 
anything  for  the  store  to  strike. 


CONCLUDING  REMARKS 

The  results  cited  in  Fig.  3-10  suggest  that  the  background  calcula¬ 
tion  is  quite  valid  in  determining  the  safe-separation  characteristic. 
Further,  even  if  the  theoretical  basis  were  unacceptable  (although  it 
seems  to  be  satisfactory) ,  the  ordinate  and  abcissa  displayed  here  seem 
tc  describe  the  data  and  separation  limits  adequately. 


SUMMARY  AND  COMPUTATIONAL  OUTLINE 

The  basic  properties  needed  to  compute  the  separation  character¬ 
istics  of  a  given  store  are  the  relative  velocity  of  a  point  with 
respect  to  the  rack  and  the  relative  acceleration  of  the  same  point  with 
respect  to  the  rack,  both  evaluated  in  rack  coordinates  at  the  instant 
of  release.  The  former  may  be  written 


A  Jk 


[*ol  ’  (Bo  +  uo  x  rp)  -  VV  -  x  ([Ro]  •  rP) 


Co. 


Data  courtesy  of  Charles  Dragowitz  of  Grumman  Aircraft  Engineering 
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A— •  FLIGHT  TESTS 


FIG.  7.  Vurtical  Separation  Characteristics  of  Mk  4  Pod  From 


A-8  FLIGHT  TESTS 
(MK-41 
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FIG.  8.  Lateral  Separation  Characterisitcs  of  Mk  4  Pod  From  A-6 


OV-IO  FLIGHT  TESTS 
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FIG.  10.  Lateral  Separation  Characterisitics  of 
Sargent-Fletcher  Fuel  Tank  for  0V-1D. 

and  the  latter,  if  Uj  -  dU/dt  at  t  -  0,  etc.,  is 


[RJ  ‘  (Ul  +  x  rp)  -  Vj’  -  %'  x  ([R0]  •  rp) 

In  the  x',  y',  z*  system,  the  relative  velocity  components  are  (Ac.Aq.At) 
and  the  acceleration  comoonents  are  (A6,A4,A2).  Because  of  the  cumber¬ 
some  nature  of  [R] ,  simplified  procedures  are  useful. 
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Vertical  Motion 

The  procedures  described  here  are  based  upon  the  assumptions: 

1.  Aircraft  motion  is  steady. 

2.  The  store  is  rigidly  connected  to  the  airframe. 

3.  The  airframe  is  rigid. 

In  this  case  the  relative  velocity  between  the  store  and  reference  plane 
is  zero  at  a  time  just  before  the  start  of  the  release  sequence.  With 
this  simplification 

A^  *  -Aoiy  '  x|cos  ip0  sin  0O  sin  0q  +  cos  0O  cos 

+Au)y  •  z  |cos  \l>0  cos  0O  sin  %'0  -  sin  0O  cos 

+  w.  (30) 

w  is  the  ejector  velocity  and  Auiy  is  the  difference  between  its  angular 
velocity  imparted  by  the  ejector  system  and  the  aircraft,  i.e.. 

Sly'  -  uiy(0).  The  value  for  z  is  -r(x),  the  radius  of  the  store  at 

position  x.  To  continue, 

A2  -  ?  |cos  <J>Q  cos  0q  +  \  pU2(0)|  j  (cZq  +  6CZ()j  -  ^r(cM0  +  SC^jjtimes 
jcos  ip0  sin  0O  sin  0q  +  cos  0O  cos  Qq 


+  9 


-  cos 


<f>0  sin  0’  +  |pu2(0)||(cXo  +6CXo)  +  ~ +  «%)!  times 

y 

jcos  \pQ  sin  9q  cds  0o  -  sin  0O  cos  9qJ 


+  6CJ  +  0  (%  +  scOlsl”  ♦. 

y 


sin  0O  cos 


(31) 
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where 


<t>0  *  aircraft  roll  angle  at  release 
6^  -  aircraft  pitch  angle  at  release 
4>’  “  aircraft  yaw  angle  at  release 


Lateral  Motion 

The  several  terms  below  have  been  defined  elsewhere,  i.e.,  Uz^  is 
the  first-order  change  in  store  vertical  velocity. 

A3  *  -  u)y  x  sin  ipQ  sin  0O  +  v  (32) 

A4  -  UX1  sin  ip0  cos  0O  +(uyi  +  wZl  xj  cos  ipQ 

+  ^UZl  -  xj  sin  ipQ  sin  0Q  (33) 

U  -  g  sin  *o  +  ^pU^(0)“  |Cyo  +  6Cy0  +  4j§  (cNo  +  «CNo)j  (34) 

y 


Axial  Motion 


A5  ■  -  Awy  x  sin  0O  -  Au>2  x  sin  \p0 


A6  «  g 


g  -  §PU2<0>§  (cXo  +  ScJ  -  ^(cKo  +  Scj]  tines 

|cos  i|)0  sin  0O  sin  0q  +  cos  0O  cos  0qJ 
' 

*  fpu2(0)|  (c2o  +  6cJ  +  +  5C„J  tines 

L  y  J 

|cos  ip0  sin  0q  cos  0O  -  cos  @o  sin  0o] 


-  sin  0q  cos  <j)(j 
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Computational  Outline 


Tc  determine  the  safe-separation  criteria  for  a  given  store,  rack, 
and  aircraft  combination  the  following  information  is  needed: 

1.  Geometric 

(a)  Store — radius  distribution,  including  fins  if  any;  length; 
location  of  store  center  of  mass. 

(b)  Rack — foot  position  with  respect  to  store  center  of  mass; 
location  of  nearest  strike  point  in  each  direction. 

(c)  Euler  angles  of  store  and  rack  at  launch. 

2.  Force  and  Maas  Properties 

(a)  Store — weight;  moment  of  inertia  about  each  axis. 

(b)  Rack — force  versus  distance  or  time,  for  ejector  foot,  if 
any,  for  each  rack  load  and  condition.  This  can  be  inte¬ 
grated  to  give  impulse.  If  the  details  of  the  rack  ejector 
impulse  are  not  known  exactly,  the  best  possible  approxima¬ 
tion  should  be  used.  Ii?  possible,  this  approximation  should 
underestimate  the  velocity  imparted  by  the  ejection  system. 

(c)  Aircraft  flight  speed  and  dynamic  pressure. 

3.  Store  Aerodynamic  Coefficients  in  Suspended  Position  at 
Aircraft  Maneuver  Conditions ,  Velocity ,  Mach  Number ,  etc. 

The  calculation  proceeds  as  follows : 

1.  Given  initial  values  of  iy,  9,  <p  and  if>',  9',  <j>',  the  Euler  angles 
of  the  store  and  reference  (aircraft)  coordinate  system,  re¬ 
spectively,  compute  the  required  angular  transformations  given 
in  Eq.  30,  31,  32,  34,  35,  and  36. 

2.  Given  rack  ejection  characteristics  and  store  inertial  char¬ 
acteristics,  compute  v,  w,  and  AiOy^  Au)z,  the  translational  and 
angular  velocity  components  of  the  store  due  to  the  ejector. 

3.  Select  critical  points  on  the  store  x  and  V(x)  and  compute 
Ap  A^,  Aj  from  Eq.  30,  32,  and  35,  and  normalize  with  respect 

to  S2g  rmax.  Ap  A^,  and  A^  are  the  required  z',  y',  and  x' 

component  velocities  of  the  critical  points. 
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4.  Next,  given  the  free  stream  aerodynamic  characteristics  of  the 
store,  develop  the  free  stream  aerodynamic  coefficients  (in 
body  coordinates)  Cz  ,  CXq,  Cy  and  Cm  ,  Cnq  for  the  aerodynamic 

conditions  («,  3*  q)  that  exist  at  t  •  0. 

5.  Develop  the  aerodynamic  interference  coefficients  (in  body 

coordinates)  6C2q,  6CyQ,  and  ^N0  ^ro®  Eq.  113-119 

or  from  wind  tunnel  data. 

6.  Compute  A2,  A4 ,  Ag  from  Eq.  31,  34,  and  36  and  normalize  with 
respect  to  g.  A2,  A4,  and  Ag  are  the  component  accelerations 
of  the  store  critical  points  in  the  reference  coordinate 
system. 

7.  Prepare  charts  of  normalized  velocity  (ordinate)  and  normalized 
acceleration  (abscissa) ,  using  the  store  radius  and  a  critical 
time  (say,  1/4  second)  to  develop  the  inclined  boundary  from 
Eq.  16. 

8.  Plot  points  previously  calculated. 

Finally,  one  should  evaluate  the  linear  "jones'  Limit"  for  stable 
store,  i.e,. 


with  the  compuced  A2 /g  >  ^2 ^Jones  * 
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SEPARATION  ANALYSIS  IN  DETAIL 


REFERENCE  COORDINATE  SYSTEM 

The  first  problem  encountered  in  studying  store  separation  is  the 
selection  of  the  most  convenient  coordinate  system.  After  considerable 
experimentation  the  following  selection  was  made.  The  origin  of  the 
reference  coordinate  system  is  on  the  rack-store  interface  waterline; 
located  over  the  store  center  of  mass ,  in  the  sense  that  it  is  perpendicu¬ 
lar  to  the  aircraft  waterline  in  the  buttline  plane.  The  plane  defined 
in  the  introduction  can  be  identified  with  the  waterline  plane  of  the 
store  rack  (note  each  store  has  its  own  plane) . 

The  x'  axis  is  forward  in  the  buttline  plane,  the  z'  axis  is  down¬ 
ward  in  the  buttline  plane  and  y*  is  normal  to  the  buttline  plane  in  the 
waterline  plane  positive  in  a  right-handed  sense.  The  velocity  of  the 
origin  with  respect  to  the  earth,  measured  in  the  coordinate  system 
described  above  is  denoted 

v*  -  vx'  V  +  Vy  j'  +  V*'  k’  (38) 

Further,  IJhe  coordinate  system  has  angular  motion  with  respect  to  the 
earth  of  ft’.  If  the  aircraft  center  of  mass  is  located  a  distance  r' 
from  the  reference  origin,  the  velocity  of  the  center  of  mass  with  respect 
to  the  earth  in  a  coordinate  system  parallel  to  the  reference  system  is 

VA/C  =  V'  +  ft’  x  r’  (39) 

Generally  V',  ft',  VA/c  are  known  functions  of  time.  These  functions 
include  the  case  of  the  maneuvering  aircraft. 


POSITION  OF  STORE 

To  define  the  motion  of  the  store,  consider  a  body-fixed  system  with 
an  origin  at  the  center  of  mass  of  the  store  and  the  axes  aligned  with  the 
principal  inertial  axes,  x  is  forward,  z  is  downward,  and  y  completes  a 
right-handed  system.  This  coordinate  system  is  located  with  respect  to 
the  earth  by  (1)  yawing  to  the  desired  heading  (+  nose  right)  (2)  pitch- 
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ing  to  the  desired  attitude  (+  nose  up)  and  rolling  (+  right  wing  down). 
These  angles  are  ip,d,<p,  respectively.  Transformation  matrix  from  the 
x  system  to  the  x'  system  is  (see  Appendix)  quite  complicated  because 
it  involves  first  transforming  from  the  store  to  earth  coordinates 
(Appendix) . 


[<J>  •  9  *  <l>]  y 


cos  i(i  cos  0  cos  i//  sin  8  sin  $  -  sin  t |i  cos  6  cos  i(i  sin  6  cos  $  +  sin  <l>  sin  $ 

»  sin  i|/  cos  8  sin  ifi  sin  6  sin  $  +  cos  iji  cos  $  sin  i p  sin  6  cos  $  -  cos  i (i  sin  (j>  (41) 


-  sin  8 


cos  8  sin  4 


cos  6  cos  $ 


and  then  to  the  prime  coordinate  system 


I<j>'  •  9'  •  ifi']--*-  •  [if>  •  0  •  d>] 


which  will  be  denoted 


Similarly,  the  linear  and  angular  velocities  of  the  store  with  respect  to 
the  earth  may  be  computed  in  the  reference  direction 

U'  -  [R]  •  U  (44) 

Thus  the  relative  velocity  between  the  store  center  of  mass  and  the  refer¬ 
ence  point  is 

V'  -  u'  -  V'  -  [R]  •  U  (45) 

This  can  be  integrated  to  give  the  relative  distance 


dr  -  dr(0)  +  yt  [V’  -  [R]  •  V]dt 


34 


NWC  TP  4995 


Note  that  In  calcula  'ni  we  have  used  the  relations 

•t. 

>o 


<Kt)  -  d>(0)  +  /  4>(t)dt 
Jo 


(47) 


and 


tan  6  sin  <}>  tan  6  cos  $ 
cos  <J>  -sin 


sin  tj) 
cos  9 


cos  4> 
cos  6 


<*x 

“y 


(48) 


A  point  P  away  from  the  origin  of  the  store  has  a  velocity  in  store 
coordinates  relative  to  the  earth 


Up  «  U  +  u  x  r, 


(49) 


Its  relative  velocity  with  respect  to  a  point  Q  in  the  reference  system 
is 


V'  +  ft'  x  rq  -  |r]  •  [u  +  w  x  rpj 


(50) 


Usually  we  are  interested  in  the  image  of  P  in  the  reference  plane, 
z'  =  0 


?'  •*  r*  = 
Q  PI 


=  [Rl  •  t 


(51) 


thus 


<¥pip  =  ~fv'  x  (tR3  *  ^p) 

oi 

-  [R]  |u  +  a)  x  rpj  J 


(52) 


By  considering  components  of  d^'p^-p  we  know  the  vertical  distance  6z' 

and  horizontal  distance  <5y'  between  the  store  and  reference  points.  The 
actual  computation  is  not  difficult,  but  is  quite  messy  because  [R]  itself 
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is  extremely  bulky.  For  our  purpose  ve  will  assume  the  aircraft  has 
pitching  notion  only,  so  i|»'  *  <f>*  -0,  and  the  store  has  rotational 
symmetry,  <p  -  0.  In  this  simple  case, 


(*) 


feet  -li  coa  8’  coa  0  +  alu  6*  aln  8 
lain  8  c.w  6 

[cca  8  coa  8  ain  8'  -  aln  8  coa  8: 


-  coa  8‘  aln  t 
c  oa  8 

-  aln  8  coa  6  aln  8' 


cos  8  sin  6 

coa  8' 

-  aln  8* 

co-  8* 

sin  8  tin  8 

cca  8  aln  8' 

1  aln  8 

*  eo*  8* 

coa  8. 

(53) 


and  so  for  a  non-ejecting  rac«>,  with  m  and  *  0 


i 

4*'  -  4*'(0)  -  j  I  Vj(e)  -  HyX  (cos  *  coa  6  coa  8'  aln  8  aln  8*  ]  ♦  fliy  sin  8  coa  8‘ 

e 

-  (uiy  +  AwyJxHcoi  v  aln  0  *in  0'  +  cca  0  cos  0*  ]  -  [Uy  +  w**)  aln  *  sin  0*  coa  0 
-  ftJz[cos  ^  sir*  0  coa  0’  -  coa  0  air.  8'J  -  [Ux  +  <Wy  +  iuiy)*  -  a*y]  tines 

(co#  it  cos  0  sir  0'  -  air  0  coa  **  1  +  w|dt  ^54^ 


4y‘  *  4y'(  ;■  -  j  |v^t)  -  [Ux  sin 


MOTION  OF  THE  STORE 

The  equation  of  store  motion  can  be  written  in  terms  of  the  body 
axis  system  in  *  '  following  form 


dlk 


m 


^  <1J2uy  +  0//u3 


FAER0  x  +  FGRAV  X 


(56) 


m 


dU 


^X^Z  +  uz^x 


-  F  +  F 

AERO  y  GRA7  y 


(57) 


ra 


dU 


—  -  UyUx  +  t^Uyl  =  F  +  F 

'  1  AERO  z  GRAY  z 


(58) 
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+  V*(I*  -  l>'>  *  MX  AERO  <59> 

+  “*“*(I*  -  W  *  «y  aero  (60) 

-  UyU^dy  -  IX>  -  M*  „  (61) 

The  gravity  forces  are 

Fx  =  -mg  sin  0 
Fy  =  mg  cos  9  sin  $ 

Tz  a  tag  cos  8  cos  <J» 

and  from  above 


(62) 

(63) 

(64) 


dOJ 

t  — 
L*  dt 


db) 


Ly  dt 


I2 


dOJ 
_ 2 

dt 


♦  »  <{>(0) 
e  =  0(0) 

'<!>  «  W) 

The  aerodynamic  forces  are  written 

rx  AERO  =  2pu2s  [  cx  +  cxz'z’  +  5CX(0)| 

Fy  AERO  =  2plj2s  ( cy  +  cyz'z'  +  ^cy(°)| 

Fz  AERO  ~  ^Pu2s  |cz  +  czz'z’  +  <$Cz(0)J 


/*t 

+  /  ('j>x  +  tan  0  sin  $  uiy  +  tan  6  cos  <p  u2)dt 
+ j  (Wy  cos  —  a>z  sin  4>)dt 
.  ft/ sin  .  cos  <j>  \  . 

+i0  “x  +  TSie  “*)dt 


(65) 

(6o) 

(67) 

(68) 

(69) 

(70) 
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where  p  *  air  density  at  release  altitude 


o  *  tau 


-i  Ei 

Ux 


Cz  ■  -C.  cos  a  -  C  sin  a 


Cx  *  C  sin  a  -  C  cos  a 


CL*Cd,Cy  ™  drag  and  side  forces  coefficient  in  unifora 

air 


zz 


*  rate  of  normal  force  due  to  nonunifora  flow 


5Cy  *  Incremental  value  of  interference  at  the  suspension 

point 


Similarly,  for  the  torques 
=  -i«n2 


M 


Mx  AERO  ‘  2P“MC£  +  Cix'2'  +  Sci<°>) 
V  AERO  >MCM  +  S,,2'  +  SCM(0>) 
*  |pU2Sd-(cN  +  Cjj  z’  +  5CN(0)j 


M 


z  AERO 


(71) 

(72) 

(73) 


These  equations  are  closely  coupled  and  quite  nonlinear.  Hence,  any 
analytic  solution  must  involve  a  simplification  of  sorts.  Consider  a 
series  solution  of  the  form 


Uz  =  Uz(0)  +  uZlt  +  uZ2  +  ... 

for  each  variable,  so 
dUz 

dr“^1  +  uz2t+ ... 

Recall  that  ultimately  we  are  interested  in  displacement,  so 

*-2 


(74) 


z  =  z(0)  +  Uz(0)t  +  Uz^  y  + 


(75) 


(76) 
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Hence*  a  first-order  solution  in  velocity  is  equivalent  to  a  second- 
order  in  displacement.  The  variation  in  the  aerodynamic  effects  with 
distance  arises  only  in  the  third-order  terms. 

The  series  for  the  several  terms  may  be  solved  to  give 

Ui.  -  Uz(0)  +  v  cos  (6*  -  e)  +  [ux(0)<Dy(0)  -  Uy(0)u>x(0)] 

+  |oU2(0)||cz(0)  +  6CZJ  +  g  cos  8(0)  cos  <j>(0)  t  +  ... 

Ux  “  Ux(0)  +  w  sin  (6'  -  8)  +  juy(0)a)z(0)  -  Uz(0)<Dy(0)j 

+  |f>U2(0)|[cx(0)  +  6cx|  -  g  sin  8(0)  t  +  ... 

Uy  *  Uy(0)  +  v  +  {[^(Oto^O)  -  Ux(0)tUy(0)J 

+  ipU2(0)j~[Cy(0)  +  6Cy]  +  g  cos  0(0)  sin  <K0)|t  +  ... 


where  w  *  vertical  ejector  velocity,  v  =  lateral  ejector  velocity. 


u*  »  o)x(0)  + 


Uiy  «  IJy(O)  +  AU)(0)  + 


0)z  =  u)z  (0)  + 


I  —  T 

“y(0)wz(0)  +  ^pu2(0)^(^(0)  +  6czj  t  +  ...  (80) 

Ao>(0)  +  Wz (°)«*c(0>  +  |pu2(0)|^(cM(0)  +  t  +  ...(81) 

-Ixi-  ly  (xix(0)My(0)  + -|pU2(0)~(cN(G)  +  5CN)  t  +  ...  (82) 

z  z 


Jx  -  Iy 


4m(0)  -  -  ,  C,  >  0. 

S 
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Au( 0)  *  angular  velocity  Imparted  by  the  ejector.  The  value  of  Cz(0), 
Cx( 0)  and  Ch(0) ,  for  example,  included  the  maneuvering  effects  of  the 
carrier  aircraft.  Thus  Cz  and  Cx  are  evaluated  at  the  aircraft  angle  of 
attack  at  the  instant  of  release  plus,  of  course,  any  Incidence  angles 
due  to  the  rack  and  the  ejector  velocity.  Furthe  C^(0)  is  not  only 
evaluated  at, the  same  angle  of  attack,  but  for  a  maneuvering  aircraft 
includes  C^6  +  Cf^a  also.  That  is. 


CM(0)  -  C^aCO)  +  CMqQ(O)  +  C^aCO) 


c(°)  -  aA/c(0)  +  i  +  »/U(0)  +  ^ 


i  *  incidence  angle  >  0. 

If  the  aircraft  were  rolling,  an  additional  term 
b'flx'/U(0) 

would  also  be  needed. 

The  Euler  angles  can  be  computed  directly  from  the  formula  given 
earlier;  for  example,  if  8  *  0, 


(J>  -  <j>(0)  +  <Ux(0)t 

+  ^“y(0)«*(0)  +^PU2(0)|^(c£  +  6C :/) 

»  - 
after  writing  cos  4>  *  cos  <+>(0)  -  sin  4> (0)  A<f>(t)  »  etc* 

6  =  0(0)  +  ju)y(0)  cos  4>(0)  +  u)z(0)  sin  <f>{0)  J t 

+  |  cos  <f(0)  ]|(i)x(0)wz(0)  J  -  |sin  4>(0)  J|o)x (0)^(0)  J 

+  |cos  4>(0)  |  - -  ojx(0)uiz(0) 


+  iu2(°)fi(cM  +  «cM) 
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A  similar  expression  may  be  written  for  tf>(t), 


*  -  HO)  + 


sin  6(0)  ,  cos  6(0)  /r.s 

cos  0(0)  uly(0)  +  co,  9(0)  “z(0/ 


t  + 


(88) 


The  second-order  terms  are  quite  bulky  and,  as  explained  below,  are  not 
needed  so  they  are  not  presented.  Again,  if  roll  is  negligible, 

6  =  <£»  b  o  and  inserting  ejection  effects 


6  =  0(0)  +  (Wy(0)  +  Awy)t  +  ipU2(0)~(CM  +  «CM)|^  +  ... 


(89) 


ip  =  ’fi(O)  +  (wz  +  A u2)  cos  0(0)  t 


-sin  0(O)Wy(O)ajz(O)  +  |pU2(0)~(CN(0)  +  6%) 


+  ... (90) 


which  is  a  great  simplification. 

With  all  these  steps  completed  we  can  compute  6y'  and  6z'  to  deter¬ 
mine  the  fate  of  the  store. 


SEPARATION  DISTANCE 

The  previous  integral  was  written  in  general  terms.  To  make  the 
calculation  more  specific  we  will  hold  Vz’,  Vy'  and  fiy'  constant  during 
the  calculation.  Hence,  we  write 


cos  0'  =  cos  Q'(0)  +  ftyt  +  ... 


cos  0*  (0)  -  6’  (0)ftyt  -  + 


(91) 


cos  0  -  cos  0(0)  -  0(O)ody(O)t 


^y2(0)  +0(0)  -  |pU2(0)f^CM  +  scMj 


+  ... 


(92) 


In  this  expansion  only  first-order  terms  will  be  retained,  which  will  lead 
to  an  expression  for  6z'  with  second-order  terms.  Thus  6z'  takes  the  form 


<5z'  =  <5z'(0)  +  Aj_t  +  A2  ~~  +  ... 


(93) 
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Evaluation  of  Aj^  and  A2  is  straightforward  but  tedious  because  of  the 
detail  contained  In  the  basic  equation.  If  we  assume  that  (a)  the  atom 
is  rigidly  connected  to  the  airframe,  and  (b)  the  airframe  is  rigid, 
then  the  relative  velocity  between  the  store  and  the  reference  plane  is 
zero  at  a  time  just  before  the  start  of  the  release  sequence.  With  this 
simplification 


cos  ip0  sin  80  sin  Q'0  +  cos  90  cos  6q j 
cos  i^0  cos  0O  sin  0^  -  sin  0O  cos  60j 


A1  =  -Aojyxj 

+AiUy2 

+  w 

where  the  subscript  o  implies  evaluated  at  time  t  ■  0.  Further, 


A,  =  ju,^  -  to y  xjjcos  ip0  sin  0Q  sin  +  cos  0Q  cos  0O| 

+  [uXi  +  tOy^zj  [cos  ipQ  cos  0Q  sin  q'q  -  sin  0O  cos  0q| 

+  [uy  +  alz1x]  jsin  ^0  sin  9o  +  9  cos  0q  cos  <j)gj 


r94) 


(95) 


Similarly,  if 


6y'  =  <5y(0)  +  A-^t  +  A^t^/2 
A3  =  -AojyX  bxn  \p0  sin  0Q  +  v 

A4  =  (uyi  +  “Zlx)  cos  ^0  +  -  00yix)  sin  <^o  sin  0O 


(96) 

(97) 

(98) 


The  several  terms  have  been  defined  above,  i.e.,  UZ1  is  the  first- 

order  change  in  store  vertical  velocity,  w  is  the  ejector  velocity,  and 
A(i3y  is  the  difference  between  its  angular  velocity  imparted  by  the  ejector 
system  aid  the  aircraft,  i.e.,  fty'-Wy(O).  The  value  for  z  can  be  r(x),  the 
radius  of  the  store  at  position  x. 

The  function  6z'(t)  represents  the  distance  between  a  point  on  the 
store  and  the  limiting  plane,  hence,  if  6z’  =  0,  the  store  is  just  at 
the  boundary  of  being  safe.  During  anv  release  or  jettison,  then,  the 
times  for  which  Sz1  =  0  are  the  critica.  period  of  the  separation  process. 
On  the  other  hand,  the  criterion  6 z’  =  C  can  be  used  to  define  those 
characteristics  of  the  store,  the  rack,  an  the  carrier  aircraft  needed 
for  safe  separation.  This  criterion  can  be  established  by  considering 
the  roots  of  <$z'  =  0,  namely  values  of  t,  defined  by  our  approximation 
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t  *  ±  y  -  2^262 '  (0) 


Here  two  criteria  suggest  themselves: 

(a)  no  real  root  for  t,  which  implies 

6z'(0)A2  >  >  0 


(100) 


(b)  no  positive  root  for  t,  which  implies 


Ax  >  0 


kj  >  6z*  (0)  •  A2  >  0 


(101) 


(102) 


To  be  able  to  apply  these  criteria,  the  expressions  for  A2  and  A2  roust  be 
put  in  simplest  possible  form.  Substituting  the  values  for  A^  and  A2  we 
find  the  criteria  for  case  a  is? 


>  [^2(0)fKo  +  <5CZJ  +  g  cos  0O 
-(^2(o^K +  {c«o]x)](co3  *0 


4>o  sin  80  sin  90  +  cos  90  cos  0OJ 


-  g  cos  4>0  +  |oU2<0)|[cXo  +  6CXqJ 

-  lpU2(0)^j^CMo  +  j^"r(x) j  ^cos  ip0  cos  B0  ain  0q  -  sin  0O  cos  0q) 

+  fpu2(0)f  (cyo  +  6Cyo)  +  ipu2(0)^  (  CNo  +  6CNo)  sin  *0  sin  d'0 

mkz  \ 


-4  w  -jr(x)Aw„  }  |  cos  cos  0  sin  0  -  sin  0  cos  0 
2  I  yoli  To  0  o  o  o 


2 

+  x|cos  i|j0  sin  60  sin  B'0  +  cos  0O  cos  0qJ 


(103) 


(ly  •  mky  -  Iz  »  rokz) 


Note  the  acceleration  ft'  x  V'  is  not  included  because  it  cannot 
act  on  the  unrestrained  store. 
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Note  we  have  used  -r(x)  for  z  because  we  are  Interested  in  the  "upper" 
surface  of  the  store.  For  case  (b)  the  first  criterion  is 

w  +  Au)y(O)  ^-r(x)  [cos  ij>0  cos  9o  sir.  ©o  -  sin  0o  cos  ©o] 

-  x[cos  ij>0  sin  90  sin  ©o  +  cos  0q  cos  ©oJ  ^  >  0  (104) 

In  examining  these  criteria,  the  number  of  variables  seems  large. 

If  one  recalls  that  ©o'-©c  *  the  incidence  angle  of  the  store,  and 
that  for  current  practice  -  0  (i.e.,  no  attempt  is  made  to  line  the 
store  with  the  average  lateral  wind),  the  number  of  variables  is  reduced. 
One  is  left  with  the  aerodynamic  forces  and  moments  acting  on  the  store 
at  the  instant  of  release,  gravity  and  the  airframe  pitch  angle.  Hence, 
for  case  (a)  one  obtains 


6 z’(0) 


>2(0)f!(czo  +  «Cz0)  cos  i  +  (cXo  +  «v)  Si-  i| 


-  j(cM0  +  ’^Mq) (x  cos  1  -  r(x)  sin  i)  I  +  g  cos  Sq 

mk 


i  w  -  Aujy(x  cos  i  +  r(x)  sin  i) 


(105) 


and  for  case  (b) 

w  -  Ac0y  (x  cos  i  +  r(x)  sin  i)  >  0 


(106) 


The  nose  corresponds  to  x  =  Zjj,  and  the  tail  x  =  -Zf. 

For  jettison  only,  i.e.,  w  and  Ao)y  =  0,  <j>  =  A/C  roll  angle  *  0, 
one  consequently  finds  for  the  (a)  case 


COS  cj)0  COS  4>0  > 


w  -  Aojy(Q)(x  cos  i  +  r(x)  sin  i) 


20(5 z'(0)  cos  <P0 


-  |pu2(0)! 


m 


(cZq  +  6cz')  "  (cMq  +  6%)  times 


d(x  cos  i  -  r(x)  sin  i) 


k2 

Ky 


(107) 
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or  for  (b) 

w  >  Atoy(0)(x  cos  i  +  r(x)  tin  i)  (108) 

For  simplicity  we  have  written 

C2»  -  Cz  cos  i  +  Cx  sin  a.  (109) 

If  the  term  in  (a)  in  braces  is  denoted  by  A,  then  if  one  plots 


and  locates  the  value  for  A,  then  A^  represents  a  safe  jettison  while 
k.2  is  unsafe. 


For  lateral  motion  we  find,  similarly. 


Note  that  the  effects  of  roll  angle  on  the  gravity  component  have 
been  included.  It  will  be  shown  below  that  the  effects  of  the  aircraft 
maneuver,  i.e.,  load  factor,  side  slip,  and  rolling  rate  are  implicitly 
included  in  the  calculation  of  the  aerodynamic  coefficient. 

AERODYNAMIC  LOADS  ON  STORE8 

The  criteria  developed  above  contain  the  aerodynamic  forces  and 
moments  explicitly.  These  applied  loads  appear  in  two  ways:  the  load 
the  store  would  experience  if  it  were  at  the  same  angle  of  attack  and 

g 

This  discussion  is  not  meant  to  be  complete,  but  rather  to  provide 
an  outline  of  some  of  the  methods  used  to  analyze  far  field  Interference. 
For  more  detail  see  Ref.  9-13. 
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side  slip  angle  (including  dynamic  effects)  in  a  uniform  field,  and  the 
load  due  to  aerodynamic  interference.  For  the  former  the  angle  of  attack 
is 


bii  ,  in' 
as  ■  aA/C  +  1  +  2U(0)  +  "u(0) 

i  *  rack  incidence  angle 

and  if  the  angle  of  yaw  is  \p0,  then  the  side  slip  angle  is 

l  ft* 

®  =  ^A/C  +  +  u(0)  (112> 

and  oi^/c  «  nW/qSC^,  n  =  load  factor.  In  steady  climbing  or  diving 

flight  n  =  cos  y0,  but  in  maneuvering  flight  where  mUooY  *  (n-cos  Yo)W  the 
value  of  n  is  not  restricted,  being  either  plus  or  minus  and  greater  than 
unity.  In  general,  this  value  for  a  is  used  to  estimate  the  aerodynamic 
load  in  a  uniform  field. 

It  will  be  assumed  that  results  of  calculation  of  uniform  field 
properties,  or  wind  tunnel  data,  are  available.  Similar  results  are 
found  for  the  side  force  and  yawing  moment  using  3  rather  than  a. 


PRIMARY  WING-BODY  INTERFERENCE 

The  aerodynamic  interference  due  to  the  slowly  varying  fields  is 


/*-£  M  C  p  •  S 

l.„2,«  n dS  ,  aMEAS  REF 

L  *  <0)J  *  “<X-!,’2)S  dx  s'(^)  -  sc-lrV 


(113) 


(+  UPWARDS) 


(k2  -  ki) 


(+  NOSE-UP) 
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■mL 


where,  assuming  an  elliptic  lift  distribution,  for  example, 


a(x,y,z) 


x  ~  -  V 

(x-xc/4)  +  z  (2tt;  AR 


(115) 


for  the  region  x  *  wing  leading  edge  to  x  =  store  nose  (a(x,y,z)  =  0 
directly  below  the  wing)  and 


ct(x,y,z) 


1  X  *  XTE 

2  2c  + 


a 


A/C 


(116) 


for  the  region  x  =  store  tail  to  x  =  wing  trailing  edge.  This  is 
essentially  Multhcpp's  procedure  (NACA  TM  1037)  for  computing  the 
induced  angle  of  attack.  If  the  flight  Mach  number  is  0.8  or  less,  the 
Prandtl-Glauert  correction  may  be  made. 


Similarly,  the  angle  of  side  slip  is 


for  the  region  from  the  trailing  edge  to  the  store  tail. 


(117) 


6  -  A 


(118) 


for  the  rest  of  the  store  length. 

Note  the  sign  on  3  is  such  that  the  local  force  is  always  toward 
the  wing  tip.  The  side  force  and  yawing  moment  are  particularly  sensitive 
to  local  configuration  changes.  The  implication  from  this  sensitivity  is 
that  in  many  cases  a  more  detailed  study  is  needed.  Such  studies  are 
available;  for  example,  Browne  and  Gallagher  (Ref.  8).  In  general,  this 
problem  has  been  studied  in  detail  and  several  references  exist  that  are 
of  detailed  use  (Ref.  5-9). 

A  second  kind  of  slowly  varying  interference  is  due  to  wing  thick¬ 
ness  t.  This  increment  is,  approximately,  directly  beneath  the  wing 

-  ( rJ 2  (119) 

and  can  be  deduced  from  the  flow  past  a  Joukowski  airframe. 
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The  unknown  Interference,  at  this  time,  is  that  due  to  the  wing  in 
the  store  flow  field  or  adjacent  stores  in  the  store  flow  field.  This 
interference  factor  Bust  be  determined  at  the  present  time,  from  wind 
tunnel  or  flight  testing.  It  may  be  optimistic  to  suggest  that  these 
complicated  interferences  can  ewer  be  computed,  but  such  a  goal  is 
surely  desirable. 
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Section  4.  CONCLUSIONS  A..I)  RECOMMENDATIONS 


The  evidence,  vLicli  is  sparse,  strongly  suggests  that  a  safe- 
separation  criteria  can  be  expressed  in  teres  of  the  relative  velocity 
of  the  store  with  respect  to  the  aircraft  and  acceleration  that  acts  on 
the  store  at  the  instart  of  release.  This  acceleration  includes  gravity, 
in  the  appropriate  direction,  the  aerodynamic  forces  divided  by  the 
weight  and  the  aerodynamic  moments  divided  by  the  inertia.  Further, 
velocity-acceleration  plane s  contain  limit  lines  determined  by  each 
individual  geometry.  On  the  vertical  plane  one  limit  line  corresponds 
to  the  requirement  that  the  store  fail  one  radius  in  a  certain  time 
(1/4  second  seems  adequate).  On  the  lateral  plane  the  line  implies  the 
store  must  fall  the  distance  necessary  to  clear  the  adjacent  store  while 
moving  laterally  the  distance  necessary  to  strike  it. 

The  proposed  criteria  have  the  advantage  that  the  information 
required  to  define  safe-separation  limits  is  the  same  as  that  required 
to  establish  air  loads  and  stresses.  Hence,  no  new  information  is 
needed  to  use  the  three  planes  (including  axial)  to  determine  safe- 
separation  requirements. 

The  proposed  velocity-acceleration  diagrams  contain  a  great  deal  of 
information  in  an  extremely  compact  way.  This  is  a  great  advantage  over 
existing  schemes.  Consider  for  example  the  flight  characteristics  of 
an  airplane  as  conventionally  represented  on  a  diagram  in  terms  of 
altitude  (h) ,  Mach  number  (ir) ,  and  load  factor,  as  shown  in  the  sketches 
that  follow — that  is,  steady  flight.  This  chart  may  be  generalized  to 
incorporate  lines  of  constant  U2/g.  This  constant  is  the  minimum  value 
associated  with  the  value  of  / /2r  ^^cr.  However,  this  value  of  &2'9 

involves  the  dive  angle,  so  an  h,M,n  diagram  must  be  constructed  for 
each  value  of  n,  and  roll  angle  <p.  A  similar  curve  must  be  constructed 
for  each  value  of  $  at  n  *  1.  Hence,  a  number  of  plots  must  be  generated 
for  the  0,  <j>  combinations,  and  this  number  is  increased  by  virtue  of  the 
number  of  different  load  factors. 

A  similar  profusion  of  charts  can  be  found  for  the  lateral  motion. 
Thus,  it  io  difficult  to  avoid  the  conclusion  that  the  dimensionless 
parameters  A]_//2rmax^,  A2 /g,  etc.,  represent  an  efficient  way  to  present 
safe-separation  boundaries . 
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On  Che  basis  of  the  results  contained  in  this  report  it  is  recom- 
oended  that  operational  methods  (such  as  a  strain  gauge  balance  in  the 
rack)  of  measuring  airloads  in  flight  be  developed.  Such  data  vould  be 
of  great  value  in  determining  detailed  correlation.  Actual  experience 
of  this  kind  in  conjunction  with  the  Aj_,  £2 »  etc.,  planes  is  also 
important  in  determining  the  sensitivity  of  score  release  characteristics 
t-  ...c  V2rlt!*",',,ns  in  store,  rack,  and  aircraft  geometry  and  per¬ 

formance.  In  thi^  -ay,  standards  for  quality  control  can  be  determined. 
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Appendix 

SOME  MATHEMATICAL  DETAILS 


GOVERNING  EQUATIONS 

The  governing  equations  are  written  in  terns  cf  a  pseudo-vector 
V  *  Vxi  +  Vvj  +  V2k  +  Ujji  +  Uyj  +  ui2k  (pseudo-vector) 


V  -  Vxi  +  Vyj  +  vzk 


(linear  velocity) 


u)  *«  +  ojyj  +  o;£k 


(angular  velocity) 


If  V(t)  denotes  the  velocity  at  t 
at  t  *  t  +  fit. 


t,  V(t  +  fit)  denotes  the  velocity 


V(t  +  fit)  -  V(t)  +  77  fit 

dt-'t 

dV 

To  calculate  we  may  use  Newton's  second  law  of  motion — 
For  translation: 


dVF 

l?  -  M  — - 
dc 


(  )  :  earth's  axes 


-*  dV  -» 
F*Mr+“xV 
at 


+  (gravity  term) 


dV  F  -  - 

—  >*  -  -  w  x  V 
dt  M 


+  (gravity  term) 
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The  aerodynamic  forces  oa  the  body  coordinate  system  are: 
Px  *  -D  cos  o  +  S  sin  8  +  L  sin  a  +  . . . ( g ) 

Fy  “  S  cos  8  +  ... ( g ) 

Tz  ■  -D  sin  a  -  L  cos  a  +  . . . ( g ) 


For  rotation: 


-»  dHE 

*-ir 


rf-^  +  SlxH 
at 


where 


^  -  W  +  tyyj  +  I2tozk  * 


(angular  momentum) 


dH_  dl  doj 

t  -»  V  -jr 

dT  -  +  ^dT  +  !^c  g 


jaJy-rf-+ j!y 


z  -*wz  dk 

+  ^CSF  +  kI^  +  ^  t 


dl  t* 

—  »  jOJj;  -  kWy 


dl 

it  -  k^v  -  iw2 


dk  -f  -f 

ae  ■  *“y  - 


Body  axes  are  principal  inertial  axes. 
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where 

Mzx  «  sin  sin  4  -  cos  sin  6  cos  4* 
MZy  »  cos  ip  sin  $  +  sin  ip  sin  8  cos  $ 
MZ2  ■  cos  0  coe  <p 
AxD  *  ^ 

AxS  “  CS  J  PS 
AxL  *  CL  \  Ps 
AyS  "  CS  \  Ps 
A2d  “  ~2  PS 

AzL  *  CL  \  PS 

where  Cp  is  drag  coefficient 

Cg  is  side  force  coefficient 
Cl  is  lift  coefficient 

and 

(  )x  denotes  X  axis  on  the  body 
(  )y  denotes  Y  axis  on  the  body 
(„)z  denotes  Z  axis  on  the  body 

lp  -  cl  i  es  \ 

Hq  -  Cb  \  pS  Al 

Nr  “  cn  \  Ps  aL 

where 

C£  ■  moment  coefficient  in  roll 
Cm  *  moment  coefficient  in  pitch 
Cn  ■  moment  coefficient  in  yaw 
Al  ■  length  of  a  store 
Bl  ■  width  of  a  store 

[Al  ■  Bl  for  a  sphere] 
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EQUATIONS  FOR  THE  EULER  ANGLES 

The  Euler  angles  0,  and  <p,  developed  In  the  conventional  order 
shown  in  Fig.  11  are  presented  here  as  functions  of  the  store .angular 
rates  u>x,  (*).,,  (uz  (Fig.  12).  The  relations  between  9,  and  in  the 
earth  coordinate  system  and  Wy,  and  in  the  store  coordinate 
system  are 

<}>  -  +  tan  G  sin  ^  toy  +  tan  0  cos  0  w2 


0  »  toy  cos  4>  —  ojz  sin  $ 


cos  4> 
cos  9 


In  matrix  form 


tan  9  sir.  4) 

cos  <j> 

sin  4> 
cos  9 


tan  6  cos 

-  sin  (j> 
cos  eh 

C08  0 


Therefore,  the  equations  for  the  Euler  angles  in  terms  oc  toy,  toy,  and  u>z 
are 


4>  -  <t>0  +  <i>  fit 

9  «*  6  +  0  <5t 
o 


'P  ■  4>0  +  5t 


where  4>c »  0O,  and  are  the  initial  conditions  for  roll,  oitch,  and 
yaw,  respectively. 
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Cl)  Yaw  +  4> 


(2)  Pitch  f  0 


*1 


(3)  Roll  -+  $ 


FIG.  11.  Euler  Angles 
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